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Abstract 

T-H Direct photons are a powerful probe to study the properties of the unique 

matter created in uhrarelativistic heavy-ion coUisions. They carry informa- 
tion on the various stages of a heavy-ion coUision. At different transverse mo- 
menta (pt), different production processes dominate the direct photon yield 
D in heavy-ion collisions. Photons at high pt can be used to study initial hard 

pLn scattering processes while photons at low and intermediate px provide direct 

I information on the hot and dense medium created in such collisions since they 

origin predominantly from jet-medium interactions and from thermal radia- 
I— I tion from the medium itself. PHENIX has measured direct-photon yields over 

^ a broad pt and energy range in different collision systems such as Au-|-Au 

^ and Cu-|-Cu, allowing systematic studies of the behavior of direct photons in 

H heavy-ion collisions. Two different methods have been used to measure the 

y direct photons. An excess of direct photons in Au-f Au collisions at low pT 

^ beyond the expectation from p+p collisions was measured for the first time via 

I I internal conversion. Furthermore, the azimuthal anisotropy parameter V2 has 

been measured. The direct photon V2 sheds light on the different processes 
contributing to the production of direct photons. The p+p data at the same 
^ energy provide a baseline for understanding the heavy-ion data, but are also 

interesting in their own right, e.g. for testing pQCD calculations or - as RHIC 
collides polarized protons - for constraining models on the gluon contribution 
to the proton spin. 

^' 

^ 1 Introduction 
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In collisions of heavy-ions or of protons, a large number of photons can be mea- 
sured. They can be roughly divided into photons from hadron decays and into 
direct photons which originate from several different sources [T]. The first source 
of direct photons is non-thermal and therefore also occurs in p+p collisions. Such 
non-thermal photons are produced in processes described by next-to-leading order 
pQCD calculations, such as initial hard scatterings or bremsstrahlung emission from 
a scattered parton. These photons dominate the overall direct photon spectrum in 
heavy-ion collisions at high transverse momenta. At lovi^ transverse momenta, ther- 
mal photons [2] are expected to be the dominant source. Such thermal photons 
are emitted by a possible quark-gluon plasma phase as well as by the hadron gas. 
Furthermore, photons can be emitted by the interaction of a parton jet with the 
medium, for example the medium can induce bremsstrahlung. Such photons are 
expected to play a dominant role at intermediate transverse momenta. Measuring 
direct photons in p+p collisions is a good test of pQCD calculations. Further- 
more, the measurement of direct photons in p+p is an important baseline for the 
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understanding of direct photon production in heavy-ion colhsions. Through their 
production in quark-gluon Compton scattering, direct photons are also sensitive 
not only to the magnitude, but also the sign of the polarized gluon distribution and 
can therefore shed light on the spin structure of protons. In heavy-ion collisions, 
photons are a good probe as they do not interact strongly and can thus traverse 
a strongly interacting medium (mostly) unaffected. As aforementioned, different 
production mechanisms are involved in creating the overall direct photon spectrum. 
The azimuthal anisotropy of direct photon production is expected to be sensitive 
to such different production mechanisms [1^, therefore it is crucial to measure both 
the nuclear modification factor and the azimuthal anisotropy of direct photons in 
heavy-ion collisions for understanding the different production mechanisms. 

2 Experimental Setup 

The PHENIX experiment is designed to study collisions of nuclei and of protons 
by looking at a large number of possible probes. The direct photon measurement 
presented in the following is performed within the PHENIX central spectrometer, 
located at both sides of the beam axis, covering tt in azimuth and r/ = ±0.35 in pseu- 
dorapidity. The detectors in the central spectrometer are arranged onion like, they 
can be subdivided roughly into three groups: tracking detectors, PID detectors, 
and calorimeters. Tracking detectors in PHENIX are the Drift Chambers (DC) , the 
Pad Chambers (PC), and the Time Expansion Chamber (TEC). These detectors 
measure the tracks of charged particles and thus they play an important role in 
measuring electrons. Particle ID detectors are a Ring-Imaging Cercnkov Detector 
(RICH), the Aerogel Cerenkov Counter (ACC), and the Time-of-Flight detector 
(TOF). The RICH separates electrons and pious below the tt Cerenkov threshold 
of 4 GeV/c. The outermost detector in the PHENIX setup is an electromagnetic 
calorimeter (EMCal) for photon and electron measurement, consisting of eight sec- 
tors, each covering 7r/8 in azimuth. Six of these sectors consist of a lead scintillator 
sandwich calorimeters, four in the west and two in the east arm. The other two sec- 
tors in the east arm consist of a lead glass Cerenkov calorimeter. Global detectors 
for centrality and reaction plane measurements include the Beam-Beam Counters 
(BBC) and the Zero-Degree Calorimeters (ZDC). A distinct reaction plane detector 
has been added to PHENIX in 2007. 

3 Measuring Direct Photons 

The large background of decay photons from hadrons such as tt^ or rj makes the 
measurement of direct photons difficult. Two different methods have been used at 
PHENIX to extract the direct photon signal from the decay photon background: the 
so-called subtraction method, and the so-called internal conversion method. The 
subtraction or statistical method has its advantages at high transverse momenta 
but fails to extract a significant direct photon yield at low px due to large sys- 
tematic uncertainties. It has been well established in different analyses |H[S]. The 
idea of this method is to subtract the decay photons from the different hadronic 
decays from all measured photons. Photon candidates are measured with the elec- 
tromagnetic calorimeter, the hadron and lepton contribution to this spectrum is 
subtracted and the spectrum is corrected for detector effects such as the geometri- 
cal acceptance and the reconstruction efficiency, leading to a spectrum of inclusive 
photons 7inciusivc- The spectrum of decay photons is obtained via simulations, tak- 
ing the measured 7r° spectrum [HIITIIH] as input, and using scaling for the other 
mesons such as 77 or w. The decay photons are then not subtracted directly but via 
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the so-called double ratio = ^, 0''°*° ■ In this ratio, common systematic uncer- 
tainties on the inclusive photon and the tt*^ spectra cancel, the direct photon yield 
can eventually be calculated as 7dir = (l - • 7inciusivo- This method has also 
been applied in Au+Au collisions to reaction-plane dependent data to obtain the 
azimuthal anisotropy parameter V2 for direct photons, and in collisions of polar- 
ized proton beams it was used for the measurement of the double longitudinal spin 
asymmetry ^ll- At low transverse momenta, the method cannot be applied as the 
remaining systematic uncertainties are larger than the small signal. Furthermore, 
at high transverse momenta, the merging of tt" decay photons has to be taken into 
account. 




Figure 1: Mass distribution of e+e" pairs for Au-|-Au minimum bias events for 
1<Pt <1.5 GeV/c. The range of the two-component fit (as explained in the text) 
is 80 < TTiee < 300 McV/c^. The black dashed curve shows /(mee) outside the fit 
range. The mentioned excess over the hadron cocktail can also be seen, comparing 
the blue fdmee) with the data points. 



The internal conversion method [S] is based on the assumption that every source 
emitting real photons also emits virtual photons that subsequently convert into low 
mass electron-positron pairs. Such internal conversion occurs for example in the tt" 
Dalitz decay, or when a source emits virtual photons directly. The mass distribution 
of such photons can be described by 

(Puee ^ 2a 1 L 4TnY ^ '^^'^^)Sdn (1) 
dniee 37r rriee \j m^^ ' 

a is the fine structure constant, S" is a process dependent form factor, S' = 1 for 

2 

direct photons in case m — > 0, or m ^ px, while S = \F{m1^)\'^{l — Tff)^ for the 

h 

Dalitz decay of tt'^ an 77. Therefore, the factor 5' = for rriee > M^, this cut-off can 

be used to separate the direct photon signal from the background. 

In a first analysis step, all e"^e^ pairs are measured in a certain px region, a cut 
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removes pairs from conversion in the detector material, while the uncorrelated back- 
ground is calculated via event mixing and the correlated background is removed 
based on like-sign pairs. The mass distribution is then corrected for the electron 
reconstruction efficiency and - in the p+p case - for the trigger efficiency. The re- 
sulting mass distribution is compared to the expectation from a simulated cocktail 
of hadron decays. An excess over this expectation can be seen in Au-f Au collisions 
as well as - though smaller - in p+p collisions. A possible source for such an excess 
are internal conversions of direct photons. In a next step, the mass distribution 
is fit with a two component function /(mee) = (1 ~ ?')/c(wee) + rfdir{mee), with 
/c(mee) being the shape of the cocktail mass distribution and fdir{fTiee) as the ex- 
pected shape of the direct photon internal conversion. Both single functions are 
normalized to the data in the low mass region rriee < 30 MeV/c^ where their shape 
is nearly identical. An example of the mass distribution together with fits for the 
different contributions is depicted in Fig. [T] The fit is done for each bin within 
several mass ranges rniow < ruee < 300 MeV/c^, r is the only fit parameter and is 
the ratio of direct photons and inclusive photons. 



4 Results 
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Figure 2: The nuclear modification factor Raa for direct photons, 7r°, and rj in 
Au-I-Au collisions at y^snn = 200 GeV. The error bars depict px uncorrelated 
errors while the error boxes show px correlated systematic errors. 




The spectra measured in p+p and Au-|-Au collisions have been extended to 
Pt — 22 GeV/c and pt = 18 GeV/c, respectively. These spectra can be used 
to calculate the so-called nuclear modification factor Raa, a measure to quantify 
possible nuclear effects. The nuclear modification factor of direct photons in Au-|-Au 
collisions at y^SNN = 200 GeV is shown in Fig. [2] together with the Raa of 7r° and 
77. In contrast to the mesons, direct photons are not suppressed at px < 15 GeV/c, 
but show an apparent suppression at highest transverse momenta. This suppression 
of direct photons is not yet fully understood. One possible explanation could be 
the isospin effect, a difference in the different cross sections c™ and ct™, leading 
to different cross sections for colliding protons or neutrons. To further study this 
isospin effect, one can look into other collision systems such as Cu-|-Gu. 
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Figure 3: The nuclear modification factor Raa for direct photons and tt" in Cu+Cu 
coUisions at ^^snn = 200 GcV. The error bars depict pr uncorrelated errors while 
the error boxes show correlated systematic errors. 
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Figure 4: The direct photon Raa as a function of pt in Au+Au collisions at ^snn = 
62.4 GeV. The error bars depict uncorrelated errors while the error boxes show 
Pt correlated systematic errors. 
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PHENIX has measured the direct photon Raa in such coUisions at ^snn = 
200 GeV for pt < 17 GeV/c, the result is shown in Fig. [s] together with the tt" 
Raa- No suppression has been observed for direct photons in such coUisions, but 
the pt range is less than in Au+Au and the uncertainties are larger at highest pt 
. As the isospin effect scales with xt = 2pt/\/s, looking at lower collision energies 
would make studies on this effect available at lower transverse momenta where the 
merging of tt*^ decay photons on the calorimeter does not yet play a role. Therefore, 
direct photons in Au+Au collisions at -^snn — 62.4 GeV have been measured as 
well, the nuclear modification factor is shown in Fig. [4] Unfortunately, the uncer- 
tainties of the measurement are too large to draw any conclusions. 




Figure 5: The direct photon invariant cross section (p+p) and invariant yield 
(Au+Au) as a function of pt- The filled points are from [5], the open points from 
[1], [5]- The three curves on the p+p data are from NLO pQCD calculations [TU] , 
the dashed lines are a modified power-law fit to p+p, scaled by Taa- The black 
curves on the Au+Au data are the Taa scaled p+p fit plus an exponential. 

At low transverse momenta {pt < 4.5 GeV/c), PHENIX has measured the ratio 
I' = 7dircct /7inciusivc in P+P ^ud Au+Au coUisious at ^snn = 200 GeV via internal 
conversion. This ratio is consistent with pQCD predictions [10 in p+p collisions 
while the experimental result exceeds a binary-scaled pQCD prediction in Au+Au 
collisions significantly. The direct photon invariant yield can then be calculated as 
^^diroct^p^-j = ^ X diV'"'^'"^'™(pT)- Fig- [5] shows the direct photon invariant cross 
section and the invariant yield for p+p and Au+Au collisions, respectively. For 
Au+Au collisions, different centrality selections are shown. The results obtained 
with the statistical method are added to the plot as well, showing a good agreement 
with the internal conversion data in the overlap region. The pQCD calculation is 
consistent with the p+p data within uncertainties for > 2 GeV/c. The same data 
can also be well described by a modified power law function {App{l + 2p'^/b)~") 
which is represented by the dashed line in the figure. The p+p curve, scaled 
by Taa, the nuclear overlap function, is significantly below the Au+Au data for 
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Pt < 2.5 GeV/c. The Au-hAu data have thus been fit with another function, an ex- 
ponential plus the Taa scaled p+p fit function {Ae'^'^^^ -\-TAA'x{^pp{^+'^PT/b)~'^)- 
The only free parameters in the fit are A and T, the inverse slope of the exponential. 
For the 20% most central events, T has been found to be T = 221 ± 23 ± 18 MeV. 
If the direct photon excess is of thermal origin, T is related to the initial tempera- 
ture Tinit of the medium, in hydrodynamical models [llj, Tinit is 1.5 to 3 times T. 
The various models differ in Ti^it and thermalization time tq, these two parameters 
appear to be related, such that Ti^it • tq ^ const. 
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Figure 6: The direct photon azimuthal anisotropy parameter V2 as a function of 
Pt for different centrahty selections. The error bars depict uncorrelated errors 
while the error boxes show correlated systematic errors. 

The direct photon V2 has been measured with the same Au+Au data set as the 
aforementioned direct photon Raa in Fig- [2| The result is shown in Fig. |6]for dif- 
ferent centrality selections. Within the uncertainties of the measurement, the data 
are consistent with a zero f 2 though the data points might slightly favor a positive 
V2 ■ But so far, the uncertainties are too large to favor or disregard theoretical 
models. 

In p+p collisions, the double longitudinal spin asymmetry ^ll has been mea- 
sured in data sets from two different PHENIX runs. The results are shown in Fig. 
[7] together with theoretical calculations for different gluon polarizations AG I12| . 
The 2006 measurement significantly improves the accuracy of the result, but still 
the uncertainties of the measurement are too large to constrain AG. 

5 Conclusion 

Direct photons have been measured by the PHENIX experiment in different colli- 
sion systems at different energies. Two different methods have been used to extract 
the direct photon signal from the large background of decay photons. Both have 
their advantages in a certain transverse momentum range. The statistical method is 
more feasible at high , while the internal conversion method allows a direct pho- 
ton measurement at low • The statistical method has been extended to measure 
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Figure 7: The direct photon longitudinal spin asymmetry ^ll in P+P collisions at 
^Snn = 200 GeV. The data from two different RHIC runs are shown together with 
theoretical expectations for different gluon polarizations AG [12) . 

asymmetries in the direct photon production, the parameter V2 has been measured 
in Au+Au collisions, and the longitudinal spin asymmetry All has been measured 
in p+p collisions. 

The nuclear modification factor Raa is consistent with unity for central Cu+Cu 
and Au+Au collisions at 200 GeV/c. In Au+Au collisions, an apparent suppression 
of the direct photon signal sets in at pt > 15 GeV/c, an observation not yet under- 
stood. The measurement of direct photons at low px shows an excess in Au+Au 
collisions above a scaled p+p measurement. An exponential fit to this excess has 
an inverse slope of T = 221 ± 23 ± 18 MeV which is consistent with hydrodynamic 
models with Ti„it w 300 - 600 MeV. 

The measurements of both V2 and All so far suffer from large uncertainties and 
cannot yet be used to support or disregard theoretical models. However, with addi- 
tional detectors and more data taken in later RHIC runs, it is likely to reduce the 
size of the uncertainties and get more significant results for those measurements. 
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